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Abstract 

A unified framework to describe the adiabatic class of ensembles in the general- 
ized statistical mechanics based on Schwammle-Tsallis two parameter (q, q') entropy 
is proposed. The generalized form of the equipartition theorem, virial theorem and 
the adiabatic theorem are derived. Each member of the class of ensembles is il- 
lustrated using the classical nonrelativistic ideal gas and we observe that the heat 
functions could be written in terms of the Lambert's VK-function in the large N 
limit. We also noticed that negative specific heat exists for certain range of values 
of the deformation parameters. In the microcanonical ensemble we study the effect 
of gravitational field on classical nonrelativistic ideal gas and a system of hard rods 
in one dimension and compute their respective internal energy and specific heat. 

PACS Number(s): 05.20. Gg, 05.70. Ce, 02.30.Gp 
Keywords: Two parameter entropy, adiabatic class of ensembles, ideal gas, hard rods, 
gravity. 



I Introduction 



A system in thermodynamic equilibrium with its surroundings can be described using 
three macroscopic variables corresponding to the thermal, mechanical and the chemical 
equilibrium. For each fixed value of these macroscopic variables (macrostate) we have 
many possible microscopic configurations (microstates). A collection of systems existing 
in the various possible microstates but characterized by the same macroscopic variables 
is called an ensemble. The thermal, mechanical, and the chemical parameters can be 
choosen between either an extensive variable or an intensive variable and so, we have, in 
total eight different ensembles. The eight ensembles are further divided into two classes 
namely the isothermal class for which the thermal equilibrium variable is the temperature 
and the adiabatic class for which the thermal equilibrium variable is the heat function. The 
isothermal class comprises of the the canonical (N,V,T), isothermal-isobaric (N,P,T), 
grandcanonical (//, V,T) and the generalized ensemble (/i,P, T). All the individual mem- 
bers of the adiabatic class of ensembles have same value for the heat function. The heat 
function is defined through the relation 

Sj = E + ^2x t X t , (1.1) 

{1} 

where E is the internal energy and the x is an intensive thermodynamic variable whose 
conjugate extensive variable is X. Defining Xi and £2 as the variables corresponding 
to the chemical and mechanical equilibrium the specific form for each ensemble, its heat 
function and the corresponding entropy are listed in the Table 1. 



Table 1: List of Adiabatic ensembles 



e 


Ensemble 


Heat function 
Sj = E + J2x e X e 

{£} 


Entropy 
S(£i,£ 2 ,J5) 





Microcanonical 
(N,V,E) 


E 

(Internal energy) 


S(N,V,E) 


1 


Isoenthalpic - isobaric 
(N,P,H) 


H = E + PV 
(Enthalpy) 


S(N, P, H) 


1 


Third adiabatic ensemble 
(M,V,L) 


L = E - fiN 
(Hill energy) 


5(m,V,L) 


2 


Fourth adiabatic ensemble 
(H,P,R) 


R = E + PV - fiN 
(Ray energy) 


5( M ,P,R) 



Of the four adiabatic ensembles, the microcanonical and the isoenthalpic-isobaric en- 
sembles described in [TjH] are well known, but the other two ensembles namely (/x, V, L) 
and (/1, P, R) introduced through References [offi] are relatively less known. But they are 
important in studying adiabatically confined systems with variable number of particles. 
For example in Ref. [7] a Monte Carlo simulation of a system of liquid palladium has 
been carried out in the (/i, P, R) ensemble. The simulation was much more convenient in 
this ensemble and the results agreed with those obtained in the grandcanonical ensemble. 
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Realizing the importance of classifying the ensembles a unified treatement of the adia- 
batic ensembles was carried out in [8J. Later on the physical realizations corresponding 
to the eight ensembles and their interrelations through Laplace-Legendre transforms was 
studied in [9]. 

Tsallis proposed a generalization of the existing Boltzmann Gibbs statistical mechan- 
ics [TU] through the introduction of g-deformed logarithm and exponential functions. The 
entropic expression in the generalized statistical mechanics was based on the g-deformed 
logarithm. Many other deformed entropies like the K-entropy [TTJ, the basic deformed en- 
tropies [12] etc., were also proposed. The equilibrium formulation of statistical mechanics 
based on the deformed entropies were developed and applied to a wide variety of systems 
like long range interacting rotators [13], relativistic gases [T^fTo'] . and systems with long 
range microscopic memory [16J. Classification of the eight ensembles in to two different 
classes namely, the isothermal and the adiabatic class, and a unified description for each 
class in the framework of generalized statistical mechanics based on Tsallis q entropy has 
been carried out in [T7] . 

Later it was shown in [IB] that an information theoretic entropy known as the Sharma- 
Taneja-Mittal entropy [T9"f2"0"] based on a two parameter logarithm was a natural general- 
ization of the K-entropy. But the functional form of the two parameter exponential does 
not exist and so the distribution could not be written explicitly. With a view to provide 
a generalization of the Tsallis q entropy, a two parameter logarithm and its inverse func- 
tion the two parameter exponential were introduced in [21]. The doubly deformed (q,q') 
logarithm and its inverse the (q, g')-exponential for a real variable x are 



exp(i-^y-« -1)J-1 



1 + ^— -,la{l + {l-q')x) 
1 — a' 



i 

l-q 



1.2) 
1.3) 



The two parameter algebra based on the doubly deformed logarithm and the exponen- 
tial was found [22] to be nondistributive in nature. Based on the definition of the two 
parameter logarithm (11. 2p the generalized entropy is 



= ^^{^{t^^ 1 -^)- 1 )- (L4) 

Under the assumption of equiprobability i.e., pi = 1/w the two parameter entropy (jl.4p 
reads: 

( om ( 1 ~ g^ ,„,i-g 



S m , = ^ ^exp [j^^ (!-5) 

Both the doubly deformed logarithm (jl.2p and the exponential f ll.3p reduce to the respec- 
tive one parameter logarithm or exponential defined in [10] . when we let either q or q 1 
tend to unity. Due to this property, the Tsallis generalized entropy can be recovered by 
allowing either of the deformation parameters to approach it limiting value. 
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A canonical formulation of generalized statistical mechanics based on the two param- 
eter entropy (II. 4p was carried out in Ref. [23] and the authors found that the distribution 
could be obtained in terms of the Lambert's ^-function. This can be formally extended 
to describe other members of the isothermal class of ensembles. In our current work 
we investigate the adiabatic class of ensembles of the generalized statistical mechanics 
based on the two parameter entropy. We provide a unified description of the adiabatic 
class and provide a two parameter generalization of the equipartition theorem, virial the- 
orem and the adiabatic theorem. All the four ensembles are illustrated using the classical 
nonrelativistic ideal gas and their corresponding heat functions and heat capacities are 
obtained analytically. For the ideal gas, the phase space volume in the microcanonical 
and the isoenthalpic-isobaric ensemble can be found easily. To compute the respective 
heat function we need to adopt the large N limit and the results are obtained in terms of 
the Lambert's W- function. A brief introduction to the H^-function is provided in Section 
II of this article. Corresponding to the (/x, V, L) and the (/i, P, R) ensembles the phase 
space volume could not be evaluated exactly and we overcome this by using the large 
N limit. Since the large N limit has been already used in the computation of the phase 
space volume no further approximations are needed in the calculation of the heat function 
and the specific heat. The effect of gravity on a classical nonrelativistic ideal gas, and 
a system of one dimensional hard rod gas has been investigated [24125] in the extensive 
Boltzmann-Gibbs statistical mechanics. The authors found that the entropy is a decreas- 
ing function of the gravitational field indicating that gravity had an ordering effect in a 
thermodynamic system. In the present work we study the effect of a gravitational field 
in the framework of two parameter (q, q') entropy. First a .D- dimensional nonrelativistic 
classical ideal gas confined in a finite region of space, subjected to an external gravita- 
tional field is investigated. The entropy and the specific heat are obtained exactly as a 
function of the internal energy. In the infinite height limit, the internal energy could be 
found only when we assume that the number of particles is very large. Next we study a 
system of one dimensional hard rods in the presence of gravity. Analogus to the classical 
ideal gas, the entropy and the specific heat are obtained as a function of internal energy. 
But in the infinite height limit the internal energy could be obtained as a function of 
temperature subject to the condition that the number of particles is very large. 

The plan of the article is as follows: Following the introduction in Section I we give 
a brief summary of the Lambert's W function in Section II. A unified description of 
the thermostatistical structure of all the four adiabatic ensembles is provided in Section 
III. In Section IV, we study all the four ensembles using classical nonrelativistic ideal 
gas, and each ensemble is treated separately in a subsection. The effect of an external 
gravitational field on thermodynamic system is studied in Section V. The first part of the 
section deals with the effect of gravity on a .D- dimensional classical ideal gas. A system of 
one dimensional hard rod gas is examined in second part of the section. The investigations 
are carried out in the microcanonical ensemble. We present our conclusions in Section V. 
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Figure 1: A plot of the Lambert's W function W(x) as a function of x. 



II Lambert's W-function: a primer 

A brief summary of the Lambert's W-function and its mathematical properties is pre- 
sented in the current section. Lambert's W-function represented by W(z) is defined as 
the multivalued inverse of the function w e w = z and satisfies the relation 

W(z) e w{z) = z, zeC. (2.1) 

The various branches corresponding to the W-function are indexed by k — 0, ±1, ±2, 
and for real z < — 1/e, the function W(z) is always complex and multivalued. For real 
z in the range —1/e < z < 0, the function W(z) comprises of two real branches namely 
Wq(z) and W-i(z) represented by the solid line and the dotted line respectively in Figure 
1. The branch Wo(z) satisfies the condition W(z) > — 1, and is generally known as the 
principal branch of the W function. Correspondingly when W(z) < — 1 we have the 
W-i(z) branch. The principal branch is analytic at z = and has the series expansion 

w-EH-^- (2 - 2) 

lb. 

n=l 

with a radius of convergence e -1 . Similarly the series corresponding to the branch W-i is 

oo 

W_ 1 (z) = Y,c n f n (z); f(z) = y/2 (ez + 1), (2.3) 

n=0 

where the expansion coefficients can be computed from the recurrence relations 

k — 1 / Cfc_2 *^fc— 2 \ C-k—\ \ ^ / A \ 

Ck = kTT {— + — ) - T ~ —V d * = 2. c i c *+w ( 2 - 4 ) 

c = -1, ci = 1, rf = 2, di = -1. (2.5) 
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The series converges for — 1/e < z < 0, which covers the whole domain of the W-\ branch. 
From (12. ip the derivative of the W function can be found 

ir < 2 > = mrww)- < 2 - 6 > 

In the current article physical requirements restricts our choice of the W function to 
the principal branch. An excellent introduction to the Lambert's VT-function and its 
applications to engineering problems is discussed in [26J. Recently in signal processing 
the role of the W-i branch was analyzed [27]. The Lambert's function has been used in 
the study of both classical [28J and quantum statistical mechanics [2Uf3"U] . It also occurs 
in the study of Fokker Planck equation in the small noise limit. 



Ill Adiabatic ensemble: Generalized formulation 

The individual members of an adiabatic ensemble have the same value of the heat func- 
tion though they can be at different temperatures. For a system in thermodynamic equi- 
librium, there are four different adiabatic ensembles namely the microcanonical ensem- 
ble (N, V,U), the isoenthalpic-isobaric ensmeble (N,P,H), the adiabatic ensemble with 
number fluctuations (/x, V, L), and the adiabatic ensemble with both number and volume 
fluctuations (/i, P, R). The microstate of a system of iV particles can be represented by 
a single point in the 2DN dimensional phase space. Corresponding to a particular value 
of the heat function which is a macrostate, we have a huge number of microstates. We 
need to compute the total number of microstates since it is a measure of the entropy. The 
points denoting the microstate of the system lie so close to each other that the surface 
area of the constant heat function curve in the phase space is regarded as a measure of 
the total number of microstates. For a system described by a Hamiltonian H the surface 
area corresponding to a constant heat function $) curve can be calculated from 

N 

X {e} ' ^ {£} 8=1 

Similarly the phase space volume enclosed by the constant heat function .fj curve is 

TV 

E(Xt, X 2 , K) = J2 I^Tdn [ [ { U + E x ^ - II d ^ dD P- ( 3 - 2 ) 
x {ey ■ JnJpi {e} i=1 

Computation of the volume of the constant heat function curve assumes significance 
because of the difficulty in calculating the area of the curve. The volume enclosed by the 
phase space curve and its surface clXCcl £1X6 related via the expression 

fi(£i,£ 2j i3) = ^-E(X U X 2 ^). (3.3) 

OX) 

Since the phase space volume is a measure of the number of the microstates of the system, 
the two parameter entropy can be directly obtained from the knowledge of Z , (Xi,X2,i5) 
via the relation 

S q , q ,(X u X 2 ,f)) = k \n q , q ,E(X u X 2 ,ty. (3.4) 
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For a given adiabatic ensemble the temperature is defined via the relation 



T _ fds^y 1 (j7(Si,s 2 ,j3)r (35) 

V J kn(X 1 ,X 2 ,S))exp^((E(X 1 ,X 2 ,S3)y-^-l)y 

Using the definition of the temperature (13. 5p , the phase space volume (13.21) and the surface 
area (13.11) . we can calculate the expression corresponding to the heat function. From the 
expression for the heat function the specific heat can be calculated through the relation 



(3.6) 



For any adiabatic ensemble the expectation value of an observable O is defined as 
1 1 C C N 

< o> = £ m™ I l°< n+ z xtXt ~ a ) n dDr - d ^ < 3 - 7 > 

V x ' ^ J > x {t} Jr * J P* {£} i=l 

The above expression is extremely useful in computing the average energy in the (N, P, H), 
(/i, V, L) and the (//, P, R) ensembles. By using the suitable Legendre transformations, 
the average energy can also be obtained from the heat functions corresponding to these 
ensembles. From the entropy the extensive thermodynamic quantities whose conjugate 
intensive variables are held fixed can be computed through the relation 

X e (X 1 ,X 2 ,Sj) = -- JLs&^fi). (3.8) 

The equipartition theorem is derived in a unified manner for all the four adiabatic ensem- 
bles in the framework of generalized statistical mechanics based on the Schwammle-Tsallis 
(q,q f ) entropy. Representing the phase space variables r» and pi (i = 1,2, ...,3N) by a 
common variable jji, we find the following expectation value 



on 

9yj 



kT (S(X 1 ,X 2 ,Sj)) 1 ~ q exp ^i_l((r(X 1 ,X 2 ,i3)) 1 -«- 1)) (3.9) 



which is a generalization of the equipartition function. We notice that the (q, q') gener- 
alized equipartition theorem is dependent on the phase space volume in contrast to the 
extensive Boltzmann-Gibbs statistics. This observation is in line with the earlier result 
obtained for the Tsallis entropy in Ref. [17] and yields the Boltzmann-Gibbs limit when 
both the nonextensive parameters q and q' are set to unity. When we let the variable 
to be the coordinate and invoke the Dirac delta function, we get a specific form of the 
equipartition theorem: 

FfH \ 

(r l p t ) = (r t F i ) = kT (E(X 1 ,X 2 ,^)) l - q 



dr 



ex p ( ^z^amu^m 1 " 1 - 1) ) ■ (3.10) 
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Similarly when we set yi to be equal to Pi the momentum variable we obtain 



ft — ) = (pii-i) = kT (E(X 1 ,X 2 ,fi)y- q exp —L((E(X 1: X 2 ,?>)) 
opi/ V 1 - Q 

(3.11) 

Through a canonical transformation, sometimes the Hamiltonian of a system can be 
written in the following form Tl = Yli^i Pf + Pi Qi- Aided by equations (I3.10p and 
(13. lip we find the expectation value of such an Hamiltonian to be 

(H) = DN kT (E(X 1} X 2 ,^)) l - q exp \l-Z£{(E(2t x , X 2 , io)) 1 ^ - 1)^ . (3.12) 

A two parametric generalization of the virial theorem could be obtained from the relation 
(I3.10p and reads: 

DN \ f\-a' \ 

Y,nPi) = —DNkT Xa,^)) 1 -* exp U^S,, * 2 , *))- - 1) . 

(3-13) 

Below we verify the adiabatic theorem for the generalized statistical mechanics based 
on the two parameter entropy. Let us consider the Hamiltonian to contain an external 
parameter a, in addition to the phase space co-ordinates. An expectation value of the 
derivative of the Hamiltonian with respect to the external parameter yields the thermo- 
dynamic conjugate variable / corresponding to the parameter a 



D v 



\ da/ mM ^wl i (f ) '<* + e - *> n * 

/• (3.14) 



/&H\ 1 x - 1 f f f&H 

da 

The description of the adiabatic ensembles given above has been illustrated through cer- 
tain examples. First we consider a classical ideal gas and provide an analytic solution to 
the specific heat in all the four ensembles by using a large N approximation. Later on we 
study a microcanonical .D- dimensional classical ideal gas under the influence of gravity 
and compute its specific heat. Finally we investigate a hard rod gas under gravity in the 
framework of microcanonical ensemble. 



IV Application - Classical ideal gas 

The Hamiltonian of a nonrelativistic classical ideal gas in D dimensions is 

i 

where p, for (i = 1,2, ...,N) represent the D- dimensional momenta of the gas molecules. 
In this section we find the phase space volume corresponding to this Hamiltonian for all 
the four ensembles. From the phase space volume we derive the relevant thermodynamic 
quantities like the entropy, the heat function and the heat capacity. 
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A. Microcanonical ensemble 



The classical nonrelativistic ideal gas described by ( 14. ip is studied in the microcanon- 
ical ensemble. In order to compute the entropy of the system, we calculate the phase 
space volume enclosed by the constant energy curve. Substituting the expression of the 
Hamiltonian in (13. 2p 



N\h DN J..L. \<-^2m 

The phase space integral is computed in the following manner: First we notice that 
the momentum integration is the volume of a DN- dimensional sphere of radius yhnE. 
Next we integrate over the position co-ordinates to obtain the phase space volume in the 
microcanonical ensemble 



V N M N 



. DN 



W,v,E) = w T{ m + 1) ^~ ( 4 - 3 ) 

where we define Ai = (2irm/h 2 ) D ^ 2 for the sake of convenience. The corresponding surface 
area enclosed by the phase space curve is 



V N M N 



. DN 



m,V,E) = — Yr ^E^-\ (4.4) 

The microcanonical entropy of a classical ideal gas obtained from the knowledge of the 
phase space volume (14. 3p reads: 



a k 



M 1 - q> 
where the factor E mc is defined as 



exp(^((~ mc £-)^-l)) 



(4.5) 



V N M N 

Emc = ah r (M + 1) • (4 - 6) 

In the limit (q,q') — > 1, we recover the extensive Boltzmann Gibbs entropy. From the 
definition of temperature (I3.5p . we arrive at the following expression 

1 _ DN u^- q ^i-aim-i ^f 1 -* 



T 



k E^E^^ exp (^((H mc E^r* - l)) . (4.7) 



Inversion of the above relation (14. 7p . to obtain the internal energy is intractable and so, 
we look into the naturally occuring large iV limit. In this limit we can safely neglect the 
factor of one in comparison with (1 — q)DN/2 and this enables us to make the following 
approximation is a valid approximation because, for very 

large values of N the deformation parameter (1 — q) should be 0(1/N) for the factor 1 
to make a resonable contribution. Using the assumption outlined above we get 



exp 



1 — Q\ 2 _•. „ „/ n „\ DN I 1 — O' ,_ _ DN , i _ \ 



1 — q J DN kT mc \l-q 
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The inversion of the above function leads to an expression of the internal energy in terms 
of the Lambert's W function as 



E 



1-q 1 fl-q' fl~q'\ 2(3 



(1-q) DN 

(4.9) 



l-q'~} mc q "\l~q "\l-qjDN 

The requirement that the entropy be concave decides the range of the deformation pa- 
rameters. From the discussion in Ref. [21 J, we notice that the entropy is concave in the 
region q + q' > 1 excluding the region < (q,q f ) < 1, where the entropy does not have 
fixed curvature sign. The argument of the W function in (14. 9p is positive in the region 
(a) for which q > 1 and q' > 1, and hence the energy in this region is characterized by 
the principal branch. Under the restriction q + q' < 1, there are two other major regions 
namely, (b) q > 1 and — oo < q' < 1 and (c) q' > 1 and — oo < q < 1 in which the argu- 
ment of the W function is negative. Though we have a choice between the Wq and the 
W-i branch, the continuity requirement on energy restricts our choice to the Wo branch. 
Thus we conclude that only the principal branch contributes in the definiton of energy. 
In all the subsequent discussions in this article, we maintain the definition of the regions 
(a), (b) and (c) as described above. The specific heat at constant volume evaluated from 
the internal energy is 

BE kd Wn(bB) i-(i- q )DN/2 

C v = — = — 01 P> aW (b8)} v-ddw , (4.10) 

V dT (1 -q)DN/2 l + Wo(b/3) 1 01 f )l ' 1 J 

where the factors a and b £1X6 cLS defined below: 

ex P 1 ~ ( 4 - n ) 



1-g'H^' ^ 1-q l \l-qjDN 
Invoking the large N limit in the final expression of the specific heat we arrive at 

r k _l 1 uu) 

V {l-q)DN/2 a(l + W (b/3))" 1 ' ) 

From the expression of the specific heat f)4.12p . we notice that it can be either +ve or 
— ve depending on the values of the deformation parameters q and q' . In the region (a) 
and region (c) the specific heat is positive, whereas in the region (b) it is negative. We 
do not recover the Boltzmann Gibbs statistics from our above calculations, because in 
the computation of energy ( 14. 9p . we have made use of the large iV limit which does not 
commute with extensive (q, q') — > 1 limit. 



B. Isoenthalpic-isobaric ensemble 

A system which exchanges internal energy and volume with its surroundings in such a 
way that its enthalpy remains constant is described by the isoenthalpic-isobaric ensemble. 
In order to calculate the thermodynamic quantities we first find the phase space volume 
enclosed by the constant enthalpy curve. The integral expression corresponding to the 
phase space volume is 

W,P,V = Emps SJ p e (^n + PV - H ) rR".dV (4-13) 

V jTl Jpi \ i / i=l 
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Using the fact that the momentum integration is nothing but geometrically equivalent 
to a DN dimensional sphere of radius \j2m(H — PV), and integrating over the position 
co-ordinates we arrive at 

w, p.h) = %t T (4 + i) £ yN(H ~ PV) (4 ' 14) 

In the next step we consider a summation over the volume eigenstates. Since the volume 
states are very closely spaced the summation is replaced by an integration. But an integral 
over the volume leads to overcounting of the eigenstates. To overcome this we employ the 
shell particle method of counting volume states developed in [52E3]- In this method we 
take into account only the minimum volume needed to confine a particular configuration. 
The minimum volume needed to confine a particular configuration is found by imposing 
a condition wherein we require atleast one particle to lie on the boundary of the system. 
All the equivalent ways of choosing a minimum volume is treated as the same volume 
eigenstate and is considered only once. Using this shell particle technique to reject the 
redundant volume states we arrive at 

/ 1 \ N u-25 n at 

Z { N,P,H) = M» (-) 3 = — + N. (4.15) 

A similar evaluation of the surface area enclosed by the curve of constant enthalpy yields 

Q(N,P,H)=M N w (4.16) 

From the knowledge of the phase space volume (I4.15P , the entropy of a classical ideal gas 
in the isoenthalpic-isobaric ensemble is found: 



s, k 



q ' q ~ 1 - q> 



exp (yt((E ieH y-« -1))-1 



(4.17) 



The Boltzmann Gibbs entropy is recovered in the limit (q, q') — > 1. The factor S ie used 
in the above equation is 

^e = M N (±Y — - - . (4.18) 



pj r(D + r 

The partial derivative of the entropy (I4.17P with respect to the enthalpy gives the tem- 
perature of the isoenthalpic-isobaric ensemble 



i = 2) kE^H^- 1 exp (^((Sfa H*r« - l)) . (4.19) 



In order to invert the above equation we assume the large N limit and this consequently 
leads to the approximation that ~ jj(i-g)s> Rewriting the equation (14.191) 

based on the assumption yields the expression 



exp 



fcr) - ^ H(1 -^ ^ (^7 (s " ■ (420) 
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The solution of the above equation yields the enthalpy in terms of the Lambert's W 
function 

1 - q' (l -q'\ p 

2) 



H 



Q 



Wn 



cxp 



l-q 



(1-9)33 



(4.21) 



1' ^ " V 1 - 1 

The requirement that the entropy should be concave, along with the fact that the enthalpy 
is a continuous function restricts our choice of the W function to the principal branch. 
The specific heat at constant pressure computed from the enthalpy (14.211) is 



C, 



dH 



k/3 W (b(3) 



(l-g)D 1 + W (bf3) 
where, the factors a and b defined below: 

1 - q 1 „ l-q' 



l-(l-q)V 

[aW (bp)} w 



a 



q' S 



q 



exp 



l-q' 

1^7 



1 



In the large iV limit the expression for the heat capacity reads: 

hp 1 



Cr 



[l-q)T> a(l + W (bp)Y 



(4.22) 



(4.23) 



(4.24) 



The specific heat is positive in the regions (a) and (c) and negative in the region (b). 
Since we adopted the large N limit in evaluating the enthalpy (14,21]) we do not recover 
the corresponding Boltzmann Gibbs result, because the extensive limit and the large iV 
limit do not commute with each other. 



C. (/i, V, L) ensemble 

The Hill energy L is the heat function corresponding to the (/i, V, L) ensemble. The phase 
space volume enclosed by the curve of constant Hill energy L is 

^ V > L) = E 7nP» / / 8 (E £l ~ "» ~ L ) II (^) 

JT * J Pi \ i / 1 = 1 

Integrating over the phase space variables and r i; we arrive at 

°° V N M N 

S ^ v ^ = J2w rWTTT (L + » N)nf - (426) 

7V=0 V 2 ^ L ) 

An exact evaluation of the summation given in (14.26]) could not be achieved. In order to 
calculate an approximate value in the large iV limit we carry out the following procedure. 
First we approximate the factor (L + fiN)~ as follows 

— 

(L + „iV)^ = (MO- (l + ^p) ' , 

, . d n (DN ( L \\ . , 

= (^AT)- eX pf— mfl + — )). (4.27) 
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For very small values of x we can make use of the approximation ln(l + x) ~ x and this 
leads to 

(L + fiN)— = (jiN) 2 exp — - . 



2 



(4.28) 



In the next step we use the Stirling's approximation for the gamma function as given 
below: 



r 



DN 



1 



DN 



cxp 



DN 



Substituting the relations (14.281) . and (I4.29P in (14.261) we arrive at 



2J(n, V, L) = exp 



D L 

~2 Jt 



N=0 



V 



N 



M 



N 



(^)T-exp(-^) 



DN 
2 



(4.29) 



(4.30) 



Now the summation over N can be carried out in (I4.3(jp . enabling us to write the approx- 
imate expression for the phase space volume 







(f 


3 exp ( 





Vli D l 2 e D ' 2 M 
(D/2) D / 2 



(4.31) 



The approximate computation used by us to calculate the phase space volume (I4.3ip can 
be considered as a first order approximation of the summation in ( I4.26|) . Similarly the 
surface area enclosed by the curve of constant L is also found 



V,L) = ^- exp 



D L 

~2 Jt 



cxp 



V [ x D l 2 e D l 2 M 



(D/2) D / 2 

The entropy of the classical ideal gas in this adiabatic ensemble found from (14.311) is 



(4.32) 



9 - k 



cxp 



w 

l-q 



-he exp 



D L 

~2 Ji 



l-q 



- 1 



- 1 



(4.33) 



where, for the sake of convenience we define: 

( V^ D / 2 e D ' 2 M 
^ e_6XP V {D/2) D I 2 

The temperature of this ensemble calculated from the defining relation (13.51) is 



(4.34) 



T 2/i he 



cxp 



D L 

~2 Ji 



l-q 



exp 



-he exp 



D L 

~2~ Jt 



l-q 



(4.35) 



Rewriting (I4.35P we immediately recognize it as the equation whose solution is given by 
the Lambert's W- function. Thus the Hill energy in terms of the W- function reads: 



2/i 



In 



l-q 1 



Wo 



l-q 



exp 



1 - q'\ 2/3/i 



Q 



D 



(4.36) 
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The principal branch of the W function is chosen based on the concavity requirement on 
the entropy and the continuity of the Hill energy. The Hill energy L is uniformly positive 
in the region (c) for which q' > 1 and — oo < q < 1, whereas in the other two regions (b) 
and (c) it is postive only when the argument of the logarithm is less than 1. From (14.361) 
the specific heat at constant volume is found 

Cv =&T = -(l-g)D/2 1 + W (b{iy b = 6XP [T^q- ) Dj2- (437) 

Subject to the condition q + q' > 1 and excluding the region < (q, q') < 1, the specific 
heat is +ve in the region q > 1 and — ve in the region q < 1. Due to the noncommutative 
nature of the extensive limit and the large N limit, we are unable to Boltzmann Gibbs 
result corresponding to the Hill energy f )4.36j) and the heat capacity (I4.37p . 



D. (/i, P, R) ensemble 

The adiabatic ensemble with both the number and volume fluctuations is illustrated using 
the classical ideal gas in this section. The phase space volume of the classical ideal gas in 
this ensemble is 



N V Jr t J Pi y i J i=1 

(4.38) 

Integrating over the phase space co-ordinates namely Pi and we arrive at 

Efr, F,R) = EEl- fTMTTT (R + fJ>N — PV)™. (4.39) 

AT=0 V V 2 ~T~ L ) 

The sum over volume eigenstates is calculated by approximating it to an integral. Since a 
direct integration will lead to overcounting of the volume states, we use the shell particle 
method of counting which was developed in [32|33] . Using this technique eliminates 
the redundant volume states and the expression for the phase space volume after due 
evaluation is 

Finally we evaluate the summation over the number of particles in (I4.40p using a large iV 
approximation. The large N limit of the (R + /xiV) 33 factor is calculated as follows 



A? 



(R + ^Nf = GuiVf (l + 7^) = (» N f ex P (® In (l + ^) ) • (4-41) 
Using the approximation ln(l +x) ~ x for small values of x we get 

(R + {iNf = (fiNf exp (^pj . (4.42) 
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Based on the Stirling's approximation the Gamma function in (I4.40p is written as 

r(D + 1) « exp(-D) . (4.43) 
Substituting (14.421) and (14.431) . we can rewrite (14.401) as follows 



E(fi, P, R) = exp 



V R 



M 



N 



Ex 
v) 



exp (£>) 



(4.44) 



where T> = D/N. Carrying out the summation in (I4.44p we get the final expression of 
the phase space volume in the large N limit 



E(fi, P, R) = exp 



PR 



(4.45) 



The computed phase space volume (I4.45P can be considered as a first order approximation 
of (14 .40 p . The surface area of the phase space curve is 



V 

tt(fi, P, R) = — exp 
H 



V R 



(4.46) 



The two parameter entropy of the classical ideal gas in the (fi, P, R) ensemble can be 
immediately obtained from the phase space volume (14.451) and reads: 



S q , q > 



q' 



exp 



l-q 



exp 



V R 



l-q 



(4.47) 



where the temperature independent factor E re of (14.471) is 



M ( fi\ v 



P \VJ 

The temperature of the (fi, P, R) ensemble is 



exp (V) 



(4.48) 



1 

T 



dR 



V 

fi 



1 — 're 



exp 



V R 

fi 



l-q 



exp 



l-q' 
l-q 



^re GXp 



V R 

/i 



1-9 



(4.49) 

Inverting the expression for the temperature in (14.491) helps us in computing the Ray 
energy, the heat function corresponding to the (fi, P, R) ensemble. The expression for the 
Ray energy of the classical ideal gas reads: 



R 



[l-q)V 



In 



i-q _ 

l-q' E 



l-q 



l-q' 

i-q 



exp 



l-q'\ 15 fi 



q J V 



(4.50) 



The continuity requirement of the Ray energy R, and the concavity conditions on the 
entropy restricts our choice of the ^-function to the principal branch. In the region (c) 
for which q' > 1 and — 00 < q < 1 the Ray energy is uniformly positive whereas in the 
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other two regions viz (a) and (b) it is postive only when the argument of the logarithm is 
less than 1. From (I4.50p the specific heat at constant pressure is found: 

OR \i kj_ 

P dT (1 -q)V l + W (b/3)' 1-q 

Analogous to the (/i, V, L) the specific heat is positive in the region q < 1 and negative 
in the region q > 1 under the condition that q + q' > 1 and excluding the region < 
(<?,<?') < 1- Further the Boltzmann Gibbs results could not be recovered from (I4.50P and 
(14.5 ip since we have used the large N limit in the computation of the density of states 
and it does not commute with the extensive (q, q') — > 1 limit. 



V Gas molecules under gravity 

In order to investigate the effect of an external field on a thermodynamic system we study 
a system of gas molecules in the presence of a gravitational field in the current section. A 
system of D-dimensional classical ideal gas under the effect of gravity is explored in the 
first part of the section. To understand an interacting system in an external field, we, in 
the later part of the section, examine a system of hard rods confined in a linear box of 
length L in the presence of an external gravitational field. 



A. D-dimensional ideal gas under gravity 

A system of N noninteracting classical point-like particles of mass m is considered in 
.D-dimensions in the presence of a uniform gravitational field g. The acceleration due 
to gravity acts along one of the dimensions labeled as the z co-ordinate and the force is 
oriented in a direction opposite to that of the co-ordinate axis. To conform to reality 
we consider that the gas molecules are located within a fixed region in space and so the 
position co-ordinates take values over a finite interval. Of the D position co-ordinates 
denoted by (i = 1,2, D), the z co-ordinate ranges over the interval (0, H) and 
the remaining D — 1 co-ordinates range over the values (0, L). The Hamiltonian of the 
D-dimensional ideal gas under gravity is 

DN 2 N 
i=l j=l 

The fact that the gas molecules are to exist within a finite region in space is indicated in 
the Hamiltonian through a potential which mimics the presence of a wall 

JO within the allowed region, 
1 oo outside the allowed region. 

Substituting the expression for the Hamiltonian (15. ip in the integral for the phase space 
volume enclosed by a constant energy curve, we arrive at 

i r r ( DN r? - \ N 

m ^ E) = wi^ / / 9 E£ + E^ + f/ -^ II d ^ d V ( 5 - 3 ) 

Jr i J Vi \i=l 3=1 / i=l 
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Integrating over the phase space co-ordinates, the volume enclosed by the constant energy 
curve is 

V N M N 



S(N,V,E) 



N\h DN T(D + 1) (mgH) N 
where &(a, b) is a finite sum defined as below: 



& D (E,mgH), 



(5.4) 



N , . 

& D (a, b) = J2(-l) k ( ) (a - kbf 0(a - 
fc=0 ^ ' 



kb). 



(5.5) 



The surface area of the constant energy curve in the phase space is 

V N M N 



Q(N, V, E) 



N\h DN T(D) (mgH) 



N 



& D (E,mgH) 



(5.6) 



where the prime over the summation (15. 5p denotes a partial derivative with respect to 
the energy. Substituting the phase space volume (15.41) in equation (13.41) we obtain the 
entropy of a D- dimensional ideal gas in the presence of gravity 



Q - _ 

1 — q 



exp 



^((S fl e D (E,mgH)) 1 -' 1 -l) ) -1 



(5.7) 



where, the factor E g is defined as 



A N 



M 



N 



D-l 



N\ h DN r (3) + 1) (mp) 



2V ■ 



(5- 



a=l 



Using the entropic expression (15. 7p and from the definition of temperature ( 13. 5 p we get 



1 -u-i-, &p(E,mgH) 
T (G D (E,mgH)r 



ex P (1—7 ®^ ( E > m 9U)f~ q - 1)) • 



(5.9) 



Looking at the structure of the relation (15. 9p given above we realize that an inversion to 
obtain the internal energy as a function of temperature is not feasible. But the specific 
heat of the system can also be computed from the knowledge of the phase space volume 
(15 ,4p and the surface area (15. 6p and the temperature (I5.9P via the relation: 



a 



V 



\9E. 



T 



g o 1 an , n_\ 



(5.10) 



Substituting the phase space volume (15.41) . area of the curve (15.61) and the expression for 
the temperature (15.91) in (15.101) the specific heat as a function of temperature reads: 



C 



v 



D k exp 



' ' {{~ 9 & D {E,mgU)) l -' 1 -l) 



l-q 

(&p (E, mgH)) q d_ 
^ 9 e' D (E,mgH) dE 



Dq(E g e D (E,mgH)) 



q-l 



]n&' D (E,mgH) - (1 - q') D 



-1 



(5.11) 



16 



From (15. lip it can be noticed that the specific heat can have both +ve and — ve values 
depending on the term in square bracket. We investigate the relevant limiting cases: (i) 
First is the g — > limit in which the gas molecules behave like an ideal gas. (ii) The 
second case is the H — > oo limit in which the (15.91) can be inverted in the large N limit 
to obtain the energy as a function of temperature. We present the relevant calculation of 
this limiting case in the discussion below: 

In the infinite height limit the limiting value of the summation (15. 5p and its derivative 

are 

lim & D (E,mgH) = E® , lim & D (E, mgH) = D E®^ 1 . (5.12) 

H— >oo H— i>oo 

Substitution of the limiting value (15. 12[) in the expression for the temperature (I5.9P leads 
to: 

^ f = Vk sJ-iJSt 1 -*)*- 1 exp (^((Sp E®) 1 -* - 1)) . (5.13) 

In the large N limit using the approximation ~ - m (Q31 . we can invert 

it to obtain the internal energy as a function of temperature in terms of the Lambert's 
W-f unction 



E 



l-q 1 (l-q' fl~q'\ B 
W exp 



q' El' 11 \l~q \l~q J © 



(1-9) » 

(5.14) 



The preconditions that the entropy should be concave, and the energy should be a con- 
tinuous function helps us to conclude that only the principal branch of the W function 
occurs in (I5.14p . The specific heat of the classical ideal gas under gravity computed from 
the internal energy is 

kB Wn(bB) i-o-m 

where the factor b was introduced in (I4.23P and a g is 

a g = ^—^-^—. (5.16) 

The large N limit of the specific heat for the ideal gas under gravity reads: 

C kP 1 (*rl7) 

y - (l-q)J) a g (l + W (bB)Y 1 • ) 

Investigating the expression for the specific heat, we find that it can be either +ve in the 
regions (a) and (c) or — ve in the region (b). Since we have already made use of the large 
iV limit we do not recover the Boltzmann Gibbs results because the thermodynamic limit 
and the extensive limit do not commute with each other. 

The one parameter limits corresponding to the temperature relation (15.91) can be 
obtained by allowing either q or q' to take the limiting value. We notice that the expression 
for the temperature is the same when we set either of the parameter to unity and reads: 

The definition of temperature (I5.18P can be inverted to obtain the internal energy and 
the specific heat in the H — > oo limit, but there is no need to invoke the large N limit. 
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B. ID hard rod gas under gravity 

The generalized statistical mechanics based on the two parameter (q, q') entropy is applied 
to a system comprising of N one dimensional hard rods under gravity Initially we consider 
a system of hard rods of mass m and length a in a finite region of the space under the 
influence of a uniform gravitational field of strenth g. The position of the centers and 
the momentum of the hard rods are denoted by the set of values (zi, z N ;p ly ...,p N ). 
Two rods with centers given by zq and z^ + \ are considered so that their edges define the 
boundaries of the region, in such way that zn+i—zq = L+a. The relevant thermodynamic 
quantities like the entropy, the temperature and the heat capacity are then computed. 
Later we assume L — > oo limit of the system and obtain the internal energy in this limiting 
case. The Hamiltonian of a one dimensional hard rod gas under gravity reads: 

N 2 N N 

n = 2m + ^ mgZi + ^2 u (\ z i _ *D + ^2 u (zi;z 0} z N+1 ). (5.19) 

i=l i=l i<j i=l 

The distance between the centers of any two particles in the system cannot be less than 
the length of the rod o and the potential corresponding to it is 

u{\ Zj - Zi \) = {r ]** (5.20) 

\Zj — Zi\ > a. 

Similarly the existence of the gas molecules within the finite region of the space is indicated 
in the Hamiltonian (I5.19P via the potential 

u= i0, z + a< Zl <z N+1 -a, 
oo, otherwise. 

Using the Hamiltonian (I5.19P in (13. 2p . the integral expression corresponding to the phase 
space volume is 

i r r ( N 2 N \ N 

^M)^ / / e lY,^ + Y, m 9Zs + U + u{\z j -z i \)-E) Y[d Zi d Pi . 

(5.22) 

To evaluate the phase space volume we first integrate over the momentum variables. In 
the next step we use the substitution = Zk — {zq + ka) for the position integrals which 
helps us in integrating over the position variables. The final expression for the phase 
space volume thus computed is 

E(N,L,E) = — — eAE -mga,mgL). (5.23) 

where the M.\ is the one dimensional equivalent of Ai. For convenience we have also 
defined the following quantities 

a = Nz + — — a, L = L — Na. (5.24) 
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The surface area of the constant energy curve is 



tt(N,L,E) 



N 



A! T (f ) (mg) 



A' 



&[ (E — mga, mgL) . 



(5.25) 



The knowledge of the phase space volume (15.231) enables us to obtain the entropy of 
the system defined through the relation (13. 4p . Combining this with the definition of 
temperature (13. 5p . we arrive at 



1 K^i-g &i(E-mga,mgL) 
T~ ^ 9h (& X (E -mga,mgL))i eXP 



where the factor E g h is defined as 



1 - q' 

1 _ ((H fl h 6i(E - mga,mgL)) l ~ q 



-.,). 



A4f 



iV 



(5.26) 



(5.27) 



The form of (I5.26P makes us realize that an exact inversion to obtain the internal energy 
is not feasible. To overcome this we use (I5.10P in conjuction with (15. 23 p . (I5.25P and (I5.26P 
to compute the specific heat as a function of the internal energy and the free length (L). 
The specific heat thus evaluated is 



3N 



2 
3 A 



k exp 



1 - q' 



(E gh 61 (E - mga, mgL)) 



l-q 



1-q 

q (E gh 61 (E - mga, mgL) ) 9_1 



(5.2* 



^gh 



1 (&i(E -mga 1 mgL)) q d , { f , 

7 =\ — 77T^ In 61 (E — mqa,mqL) 

e; (E - mga, mgL) dE 1 v y ' y ' 



3N 



Analyzing (I5.28P we realize that both positive and negative values of specific heat are 
permissible depending on the argument within the square bracket. There are two pertinent 
limiting cases namely (i) the g — > limit wherein it becomes a system of hard rods moving 
in a length L, and (ii) the L — > 00 limit. In the second limit we present the expression 
for temperature below: 



3N 



= — kE lh 9 ( E - m 9®) {1 ~ 



T 



<!)'- 



exp 



l-q' 



z. gh (E - mga) 2 



l-q 



(5.29) 

which can be inverted in the large N limit by assuming (3 A/2) — 1 w (3 A/2). This 
approximation yields the internal energy in terms of the Lambert's jy-function 



E 



l-q 1 



1 - q' S 1 " 9 



Wn 



1 - g' 
1^7 



exp 



l-q 



2/A 
3Ay 



(1-9) 3N 



(5.30) 



The stipulations that the entropy should be concave and the energy should be a continuous 
function limits us to the principal branch of the VT-function. The specific heat of the hard 
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rod gas in the large N limit can be obtained from (I5.30p and has the same regimes of 
positive and negative specific heat values as the classical ideal gas. 

The one parameter limits of the temperature (15.261) is investigated. It can be noticed 
that whether we set q or q' to unity is not relevant and we always arrive at 

1 _ i- 9 6[(E -mga,mgL) 

T ~ 9h (6i(£ - mga,mgl)Y' 1 ' ' 

the temperature relation corresponding to the Tsallis g-entropy. Though the internal 
energy and the specific heat can be obtained from (15.311) in the L — > oo limit, there is no 
need to use the large N limit unlike the two parameter case. 



VI Conclusions 

In the current work we investigate the adiabatic class of ensembles in the framework of 
generalized statistical mechanics based on Schwammle-Tsallis (q, q') entropy. We do not 
study the isothermal class of ensembles since the canonical treatement carried out in [23] 
can be extended to the other members of the class which includes the isothermal-isobaric, 
the grandcanonical and the generalized ensembles. We provide a unified description of 
the adiabatic class of ensembles which includes the microcanonical (N, V, E), isoenthalpic- 
isobaric (N, P, H), the (/i, V, L) ensemble, and the (/i, P, R) ensemble. A generalized form 
of the equipartition theorem, the virial theorem and the adiabatic theorem are obtained. 
We investigate the nonrelativistic classical ideal gas in all the four ensembles. In the 
microcanonical, and the isoenthalpic-isobaric ensemble the entropy could be found for an 
arbitrary number of particles. Using the large iV-limit, the heat functions obtained in 
terms of the temperature are expressed through the Lambert's I^-function. From the 
heat functions the respective specific heats are evaluated. To the best of our knowledge 
an exact evaluation of the phase space volume in the (/i, V, L) and the (//, P, R) ensembles 
has not been done so far. We in our current work assume a large N limit and compute 
the approximate phase space volume. From the phase space volume the entropy, the 
heat function and the specific heat of the classical ideal gas is found. The heat function 
and the specific heat are obtained in terms of the Lambert's P^-function without any 
further approximation. The preconditions that the entropy should be concave and the 
heat function should be a continuous function of the deformation parameter restricts our 
choice of the VT-function to the principal branch. The two parameter entropy is concave 
only in the region where q + q' > 1 excluding the zone < (q, q') < 1 where the entropy 
is not always concave and we analyze the specific heats only in the above mentioned 
region. For the microcanonical and the isoenthalpic-isobaric ensemble the specific heat 
is positive in the region for which both q and q' are greater than one and in the region 
where q' > 1 and — oo < q < 1. But it is negative in the region q > 1 and — oo < q' < 1. 
In the regime where the entropy is concave, the heat capacities of the (/x, V, L) and the 
(/!, P, R) ensembles are +ve when q > 1 and — ve if q < 1. The microcanonical specific 
heat of classical ideal gas and a system of hard rod gas confined in a finite region of space 
and subjected to gravity is also analyzed. The entropy and the specific heat are initially 
calculated as a function of internal energy exactly. But when we assume the height to 
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be infinite and the number of particles to be very large, the internal energy could be 
obtained as a function of temperature in terms of the Lambert's VT-function. Again only 
the principal branch of the ly-function contributes to the internal energy and the specific 
heat. In the case of the hard rod gas, the free length L and the factor E — mga plays a 
role analogus to the height and the energy in the classical ideal gas. 

The investigation of the adiabatic ensembles by using the large N limit developed in 
the current article may be extended to other entropies like the k entropy [UJ and the 
two parameter Sharma-Mittal-Taneja entropy [T9|20] . This may help us in knowing the 
similarities and differences in the thermostatistical structure of different entropies. Such 
an understanding may help us to choose the generalized entropy which may describe a 
given complex system in the best possible manner. Further, construction of a Laplace 
transform based on the two parameter exponential may be of great help in establishing 
the connection between the isothermal and the adiabatic class of ensembles. But this is 
beyond the scope of the current work and will be pursued in future. 
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